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ON THE COMPUTATION OF THE ECCENTRIC ANOMALY 
FROM THE MEAN ANOMALY OF A PLANET. 



BY DR. J. MOERISON, NAUTICAL ALMANAC OFFICE, WASH., D. C. 

The relation between the eccentric and mean anomalies is given by the 
equation 

m = E — e sin ^, (1) 

in which m is the mean anomaly, E the eccentric anomaly and e the eccen- 
tricity of the orbit. 

In computing a series of values of E, it is not necessary to go beyond 
180°, for any mean anomaly 360° — m corresponds to the eccentric anomaly 
360° — E. Since E is greater than m, let us put 

E = m -\- X, 
where x is generally a small angle depending of course on the value of e. 
In the case of the orbit of Mercury, whose eccentricity is .2056, the max- 
imum value of a; is less than 12°. Substituting in (1) we have 

m = m + X — e sin (m -f- •'*') 
or X = esin(ni + ^) 

= e sin m cos x -}- e cos m sin x 

= esinm^l— ^ + . . .^ + ecosm \x — ,^^ + ...j 

Neglecting, for the first approximation, the third and all higher powers 

of a;, we have 

1 . 2/1 — ecosmX o a /o\ 

x^+ -I ^- la;— 2r=0, (2) 

e\ sinm / 

whence x = v^2.tanj^, (3) 

where tan ^ = e t/2 . , ^'""^ , (4) 

1 — e cos m 

in which |/2 is to be taken with the positive sign and d always less than 
90°. 

The computation of the second member of the last equation will be much 
facilitated by the use of Zech's addition and subtraction Tables. Hence 
we have E = m-}-x which will never diifer more than ±15" from the cor- 
rect value of E, even when the eccentricity is as large as that of the orbit of 
Mercury. Let us represent this first approximate value of E, viz., m -\- x, 

by £?i, that is, let 

^j = m + X. (5) 

Substituting E^ for jB in (1) we have 
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m-^ = Ej^ — esin^i, (6) 

which subtracted from (1) gives 

m — m-i = E — JS'i — e (sin E — sin Ey) 

= E—Ei—2e cos 1{E + E^) sin ^{E — E^) 
= E— Ej_-~e{E—E^) cos H^ + -^i) 
= {E — -Ei)(l — ecos^i) very nearly; 



whence we have 



^ 1 — e cos J?i' ^ ' 



and E = E-i -{- i=-, for the second approxima- 

^ 1 — ecos^i' ^^ 

on. Therefore we have finally 

E = m + X + 1=^. (8) 

1 — e cos ^1 

If we require a third approximation, we may repeat the last operation 
with the corrected value of E as given by (8), but it will seldom or never 
be necessary as the last equation will generally give E within 0".01. 

The second approximation given by (7) is substantially the same as that 
given by Gauss in his Theoria Motus, Art. 11, but the method here em- 
ployed for obtaining it is much easier and more direct. 

The second correction (7) is easily computed by the aid of Zech's Tables 
before referred to. 

We will now test our formulse by the following example : 

Given m = 143° and e = .2056, find E. 

logi/2 = 0.1505150 
log e = 1.3130231 = l'3130231 

sin m = 9.7794630 cos m = 9.9023486 n 

Co. log (1— e cos m) = 9.9340098 (Zech) 9.2153717 n 

tan e = 9.1770109 (? = 8° 32' 54".8 

tan|(? = 8.8735468 

l/2 = 0.1505150 

cosecl" = 5.3144251 



logic" = 4.3384869 x = 6° 3'21".525 

Ej^ = 149° 3'21".525 
esinE, = 6° 3'26".23 

mj == 142°59'55".295 
m — m, = 4".705. 
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log (m — mi) = 0.6725596 
log(l — ecos El) = 0.0705317 (Zech). 

log 3".9997 = 0.6020279. 
.•.E= m + X + 3".9997 

= 143° + 6° 3'21".525 + 3".9997 

= 149° 3' 25".5247. 

Check. e sin E = 6° 3' 25".52 

.•. m = E — <? sin E 

= 143° 0' 0".0047. 
[f the second correction be taken at 4", which it is very nearly, E will 
be found to be 149°3'25".52 exactly. 

In computing a series of values of E, the labor may be lessened a little 
by preparing the constant logarithms, viz., log e\/2, log ]/2 cosec 1" and 
logecosec 1". 



PLANE TRiaONOMETR Y BY QVA TERNI0N8. 



BY PHOF. DE VOLSON WOOD, HOBOKEN, N. J. 

A QUATERNION may be expressed under a variety of forms. Thus if a 
and h are the tensors respectively of the unit vectors a and /9 we have (see 
Analyst, Vol. VII, pp. 124 and 127) 

q = \i' (1) 

= Tq.Uq (2) 

= 8q+Vq (3) 

= Tq [SUq + VUq). (4) 

Each of these forms has special advantages for the solution of certain 

problems, but in this article we will make use of the third one. 

Two quaternions are equal when the elements of one equal I'espectively 
those of the other. Thus, if 

q = q', 
we have 

Sq + Vq = 8q' + Vq', (6) 

and hence from the definition, 

8q = Sq', and Vq = Vq'. (6) 

The principle in algebra corresponding to this is that where an equation 

is composed partly of real and partly of imaginary terms. Such an equation 



